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Abstract 


We investigate the linear stability of static coronal-loop models 
undergoing thermal perturbations. The effect of conditions at the loop 
base on the stability properties of the models is considered in detail. We 
consider the question of appropriate boundary conditions at the loop base 
and conclude that the most physical assumptions are that the temperature 
and density (or pressure) perturbations vanish there. However, if the base 
is taken to be sufficiently deep in the chromosphere, either several 
chromospheric scale heights or several coronal loop lengths in depth, then 
the effect of the boundary conditions on loop stability becomes negligible 
so that all physically acceptable conditions are equally appropriate. For 
example, one could as well assume that the velocity vanishes at the base. 

We calculate the growth rates and eigenmodes of static models in which 
gravity is neglected and in which the coronal heating is a relatively 
simple function, either constant per-unit mass or per-unit volume. We find 
that all such models are unstabl e with a growth rate of the order of the 
coronal cooling time. The physical implications of these results for the 
solar corona and transition region are discussed. 


I. Introduction 


For some time now static models of coronal loops (e.g. , Rosner, 
Tucker, and Vaiana 1 978 j Craig, McClymont, and Underwood 1978; Vesecky, 
Antioohos, and Underwood 1979) have been widely used to interpret 
observations of both solar and stellar coronae (e.g. , Bonnett and Dupree 
1980, Orrall 1981). The key assumption in these models is that the energy 
input to the corona is constant in time so that a static solution to the 
relevant equations is possible. However, even if a static solution is 
mathematically possible, it will not be physically realizable unless this 
solution is also stable to small amplitude perturbations. It is well known 
that the coronal plasma is susceptible to a thermal instability due to the 
form of the dependence of radiative loss rate on temperature (Field 1965). 
Therefore, in order to assess their degree of validity, it is necessary to 
determine whether the static models are thermally stable or unstable. 

The linear theory for the thermal stability of static coronal- loop 
models has been investigated by a number of authors, but with differing 
conclusions. In his original work, Antiochos (1979) concluded that the 
static models were thermally unstable. Similar results were obtained by 
Hood and Priest (1980). However, Chiuderl, Einaudi, and 
Torrioelli-C'.i amponi (1981); Craig and McClymont (1981); and McClymont and 
Craig (1981a, b, c) have found that the models are either stable or that the 
growth rates for instability are too small to be physically significant. 
The origin of this difference in the results of the two sets of authors is 
in their treatment of the base of the loop models. Antiochos (1979) and 
Hood and Priest (1980) have not included cool material, T X 10 5 K, in their 
model, for which the form of radiative loss curve (e.g., Raymond, Cox, and 
Smith 1976) favors linear stability (Field 1965). In addition, Antiochos 
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has considered only perturbations with a vanishing first-order heat flux at 
the base. Chiuderi, Elnaudi and Torrieelli-Ciamponi} and Craig and 
MoClymont argue that the growth rates for instability are very sensitive to 
these assumptions so that the models can be effectively stabilized either 
by including chromospheric material in the model (Craig and MoClymont 1981) 
or by changing the boundary conditions so that the temperature 
perturbation, Instead of the heat-flux perturbation, is assumed to vanish 
at the base (Chiuderi, Elnaudi, and Torrieelli-Ciamponi 1981). 

The reason for this apparently high sensitivity of the models to the 
base conditions lies in one of the well-known properties of the static 
models: the fact that the magnitude of the conduction and radiation terms 
in the energy balance are approximately equal throughout the loop (Vesecky, 
Antiochos, and Underwood 1978). As discussed by Field (1965), we expect 
thermal instability in coronal plasma whose structure is such that 
radiation dominates. On the other hand, we expect stability if the 
structure is conduction dominated, slnoe conduction always acts to damp out 
any temperature perturbation. Hence, the static models, whose structure is 
such that these two terms are comparable, are perched near the boundary 
separating stability from instability. Even seemingly minor changes on the 
form of the base or of the boundary conditions at the base can be 
sufficient to push the models to one side or the other of this boundary. 

It is clear, therefore, that the proper treatment of the base is 
critical for determining to what extent the models are either stable or 
unstable. In the next section we discuss the question of the base 

conditions in detail. In Sections III and IV we calculate growth rates and 
elgenmodes for various static models and for various boundary conditions. 
In the final section, V, we discuss the implications of our results. 




ii. Boundary. Conditions 


As an aid to determining the appropriate boundary conditions, let us 
first write down the perturbation equations for a one-dimensional loop 
model. Assuming that all physical variables are of the form: 

f(s,t) = f q( s ) + e Vt f^s) , (1) 
and assuming that the growth or damping time 1/v is slow compared to the 
sound travel time so that the acceleration terms can be neglected ip the 
momentum equation, we obtain (Field 1965, Antiochos 1979) J 
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A ds 
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and 




» 


(5) 


where: A(s) is the cross-sectional area of the loop; gjj (s) is the 

component of gravity parallel to the loop; F-^ is the perturbed heat flux, 
which for the Spitzer (1962) conductivity is given explicitly by: 


F, = - 10 


-6 d_ 
ds 


(t q 5/2 T,) 


( 6 ) 


constitutes the perturbed energy sources and sinks of the plasma. For 
optically thin radiative losses and for a coronal heating function that 
depends only on density and temperature or, equivalently, pressure and 
temperature, -C^ is given explicitly by (Antiochos 1 979 ) s 
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(7) 


where A (T) Is the radiative loss coefficient (e.g. , Cox and Tucker 1969) 
and e(p,T) la the coronal heating rate. The equilibrium profiles are, of 
course, given by (e.g., Veseoky, Antiochos, and Underwood 1979)*. 


d_ 

ds 
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A(T Q ) + e(P Q , 


T 0> 


(8a) 

(8b) 


Note that we have neglected in equations (2) - (8) effects due to 
changes in ionization equilibrium or to optically thick radiation, both of 
which are clearly important in the chromosphere. Hence, this model cannot 
be used to determine the stability of the true solar chromosphere. We use 
this model chromosphere only for investigating the effect of having 
thermally stable material at the base of the corona and transition region. 
Since our results will turn out to be insensitive to the detailed structure 
of this region, we expect that our analysis will also be applicable to the 
solar chromosphere as long as it is thermally stable as well. 

It is evident from the equations above that there are four independent 
spatial derivatives in the problem, which Implies that four boundary 
conditions are required. This is to be expected since the full nonlinear 
equations also require four spatial boundary conditions for a unique 
solution (e.g., Richtmyer and Morton 1967)* The usual situation is that at 
end of the loop two boundary conditions are specified. One condition 
generally relates to the thermal properties of the base; for example, with 
respect to the perturbation the loop base may act as a thermal bath, in 
which case the appropriate condition is that the temperature perturbation 
vanishes there, T^ = 0; or it may act as a thermal insulator so that the 


The other condition 


perturbation heat flux vanishes there, & 0. 

generally relates to the inertial properties of the base; for example, it 
may act as a rigid wall with respect to the perturbation so that the 
velocity vanishes, v lb B °» or a "free sr face," so that the pressure 

perturbation, P s 0. Most of the discussion in the literature has 

lb 

concerned these four conditions, and the dlsoussion in this paper will also 
concentrate on these four. Let us emphasize, however, that there are an 
infinite number of other possible boundary conditions and that, lacking 
some physical justification, there is no reason to single out the four 
above. ^ priori , all conditions are equally acceptable. 

There has also been a considerable discussion in the literature on the 
distinction between symmetric and antisymmetric modes. However, this 
distinction is spurious. It is due only to the special class of models 
considered by the previous authors (Antioohos 1979} Hood and Priest 1980; 
Chluderi, Einaudi, and Torricelli-Ciamponi 1981; Craig and McClymont 1981). 
In their models gravity is neglected, and the loop area and the coronal 
heating do not have any spatial variation. Under these simplifications, 
the static equations are autonomous, and the equilibrium models are 
symmetric about the loop apex. Hence, the solutions to the first-order 
equations, i.e. , the normal modes of the loop, are either purely symmetric 
or antisymmetric about the apex. These results will clearly not hold for a 
realistic equilibrium model that has no special symmetry properties. For 
such a model each choice of four boundary conditions will, in general, 
determine a complete and distinct set of eigenvalues and eigenfunctions. 
Within each set all the modes are equally valid. Hence, in this paper we 
will avoid any discussion in terms of symmetric or antisymmetric modes; and 
since the most important mode from the viewpoint of thermal instability is 
the lowest order one, we will concentrate on the fundamental. 


The question remains as to what choice of boundary conditions best 
represents the physical situation at the base of a solar loop. Ideally, 
one would like conditions at the loop base not to affect the stability 
properties of the corona and transition region. The standard procedure for 
achieving this is to place the base sufficiently far that over the time 
scale of any instability no effects can propagate to the base. However, 
the sound travel time In the corona, chromosphere and even photosphere is 
very rapid compared to typical thermal time scales, such as the coronal 
cooling time. On this time scale, material at the loop base is always in 
contact with the coronal plasma, even if the base is assumed to be in the 
photosphere. Note that our neglect of the acceleration terms in the force 
equation is tantamount to assuming an infinite sound speed and, 
consequently, a negligible propagation time scale throughout the loop. For 
such a model it is not obvious, at least to us, whether the position of the 
base can ever be chosen so that conditions there have an insignificant 
effect on the behavior of the corona. We show below that, fortunately, 
even with the rapid sound speed, the base has a negligible influence on the 
ooronal stability if it is placed sufficiently deep. 

In order to determine the physically valid boundary conditions to 
impose at the loop base, we need a definite model for this base region. As 
a simple model let us assume that the chromosphere below a coronal loop is 
isothermal at temperature T b , is in hydrostatic equilibrium, and is 
strictly stable. In addition, we assume there is no real end point to a 
solar loop, so this model base region extends below the corona and 
transition region to arbitrarily large depths, i.e, , either to many 
chromospheric scale heights or to depths large compared to the coronal loop 
length. Finally, let us neglect the curvature of the loop for this 
analysis and set gii (s) to be a constant. 


■I " 






The form of the boundary conditions follows simply from the 
requirement that, since we are interested in coronal disturbances, the 
perturbations should vanish or at least remain finite in the chromosphere. 
Combining equations (3) and (5), an expression can be obtained for the 
pressure perturbation P in terms of the temperature perturbation T : 


¥ s) ¥°> /\_ds_ 

l 0 (s) V 0) " J T 0 H( V 


( 9 ) 


where H (T Q ) is the gravitational scale height at temperature T^: 


H(T) s 


2kT 



( 10 ) 


and s = 0 is some reference point in the loop (the top is convenient). We 
are defining s to increase downwards. It is evident from (9) that in order 
for (s) to remain finite with increasing depth in the chromosphere, i.e. , 
s -> «, the temperature perturbation must vanish with depth -> 0. In 
faot, must decay quite rapidly with depth, on the size scale of the 
chromospheric scale height H(T^) or leas, because the equilibrium pressure 
P 0 grows exponentially on this size scale. Note that we are excluding the 
possibility of modes that oscillate indefinitely with depth. We are 
interested primarily in the lowest order mode; and as we will see below, 
this usually has no zero crossings. 

We conclude, therefore, that the proper thermal boundary condition is 
that 0. This result is to be expected physically. It follows from 
the fact that any perturbation of the temperature also implies a 
perturbation of the gravitational soale height. If the temperature 
disturbance extends more than a scale height or so deep into the 
atmosphere, a very large pressure disturbance must result due to the change 
in scale height. Hence, the only physically acceptable temperature 
perturbations are those that vanish on the order of a scale height or less. 


Now, to determine the Inertial boundary condition, we integrate the 
continuity equation (2) to obtain an expression for the partiole flux n^V^ 
in terms of the density perturbation n^: 

s 

Oq(s)V^(s) = UgCO) V 1 (0) -v/ ds (11) 

As in the situation above, it is evident that in order for n^V^ to remain 
finite as s ■>“, n^ must vanish with depth. Again, this is to be expected 
physically. If the density disturbance extends over a large depth, then 
the change in the total number of particles in this region roust be large, 
and consequently a large flux of particles into or out of the region must 
result. It is evident from (11) that in the chromosphere n^ must vanish 

faster than L /a, where L. is the size scale of the coronal portion of the 

c c 

loop, i.e. , the region where we exDect a finite n^. Note that the density 
perturbation may extend to much greater depths than the temperature one. 
We found above that Tn must vanish at least as fast as exp(-s/H b ), where 
is the gravitational scale height in the chromosphere. For typical solar 
loops 1^ « L c . 

We conclude that the proper boundary conditions at the loop base are 
that T lb = 0 and n lb s 0. At any finite depth this is equivalent to the 
set T lb s 0 and P lb = 0. Although the temperature, density, and pressure 
perturbations vanish, the mass flux does not. Using (2) and (3)t the 
perturbed mass flux can be related to the pressure perturbation: 

"o v i ■ n 0 <0) V°> - Sijj - V 0)) • < 12) 

Since the pressure perturbation vanishes with depth, the mass flux will in 
general tend to some finite constant, corresponding to a steady-state flow 
in the chromosphere and below. This is to be expected on physical grounds 
as well. We do not expect the temperature and density structure deep in 


the chromosphere or photosphere to be affected by a disturbance In the 
corona. However* If as In our problem* the time scale of the disturbance 
is long oompared to the sound travol time* then a steady flow must occur at 
large depths to support art}; coronal flows. Of course, the velooity itself* 
Vj , vanishes with depth since the density increases; however, the mass 
flux stays constant. This is exactly the situation that we found for 
the case of evaporative cooling of coronal loops (Antioohos and Sturrock 
1976). Thore again* the coronal evolution oocurred on time scales long 
oompared to the sound travel time* and it was found that a steady flow was 
set up at the loop base due to evaporative motions in the corona. This is 
also the situation for the photosphere at the base of an open magnetic 
field region such as in a coronal hole* Is these regions coronal heating 
results in a steady solar wind flow. ^L’he presence of this flow has no 
observable effect on photospherlc temperatures and pressures* but it 
clearly implies the existence of a mass flux in the photosphere equal to 
the coronal mass flux. 

Let us now examine the effect of using different boundary conditions 
at the loop base, in particular* the rigid-wall one. We wish to find out 
how the coronal perturbation depends on the base conditions* for example, 
how P 1 (0) depends on n ob v lb * From (12) it is evident that specifying the 
base mass flux is equivalent to specifying the base pressure perturbation* 
but the relation between P.^0) and P lb is given by (9)i 


V°> ■ 


y°> 

P 0b 


p ib + V 0) 


flL 
J T o 


ds 

Hey 


(13) 


0 

Equation (13) is the key result. The effeot on the coronal perturbation 
P^O) of assuming a non- vanishing P^ b is given by the first term on the 
right-hand side of (13). This term is of order P n (0)/P n , for a temperature 
perturbation that vanishes rapidly in the chromosphere so that the pressure 


A ■ ■ 


.* i ■ ft* ♦ 4» ► » 


perturbation there is of order the coronal one (i.e. , (o ) does not 

diverge oxponentially in tho chromosphere). For a base that is several 
scale heights deep, P Q <0)/P 0b is negligible; therefore, the form of the 
coronal perturbation is insensitive to the boundary condition on v> 
equivalently, V^. 

The situation is clarified by Figure 1, where for a given T^s) we 

plot the resulting form of P-^s) and n Q V 1 (s) for both rigid-wall V lb s 0 

and free-surface P„ , s 0 boundary conditions. Since this is for 

lb 

illustrative purposes only, we have taken an extremely simple equilibrium 
loop model in which we neglect the effects of gravity in the corona and 
assume that the transition region is vanishingly thin. The corona is, 
therefore, Isothermal and isobaric and is defined to extend from s = 0 to 

s b L . Below this corona is an Isothermal chromosphere in hydrostatic 

c 

equilibrium and with a temperature 100 times lower. The base of the loop 
model is taken to occur 10 chromospheric scale heights deep, 
L » L c i 10 H b . For simplicity, the temperature perturbation, Tj/Tq* is 
assumed to be constant in the corona and to decay exponentially in the 
chromosphere on a size scale of H^/IO. 

The forms for P^s) and n^V^(s) obtained from (9) and (11) are plotted 
in Figure la and 1b, respectively. We have taken the velocity to vanish at 
the top, V^(0) = 0, which would be the case for the lowest order mode of a 
symmetric loop model. On the scale of this figure, the modes with s 0 
and those with V^s 0 are indistinguishable except within a few scale 
heights of the base. In this boundary region the rigid wall mass flux 
exhibits a sudden deorease to zero, and P-j^ exhibits a sudden increase back 
up to its coronal value, i.e., P lb = P-^O). Note that the magnitude of 
this jump does not change with position of the base. No matter how deep 
the rigid wall is placed, there will occur a small boundary layer at the 


wall where the pressure jumps back up to its coronal value. The 
physical reason for this is straightforward. Flows in the corona will 
cause material either to pile up or to evacuate from the wall. Hence, 
there must exist a pressure gradient near the wall of the same order as the 
coronal gradient in order to decelerate or accelerate the flows. However, 
over most of the chromosphere the flows are in a steady state, and the 
pressure perturbation is negligible. 

It is clear from these results that the appropriate boundary condition 
to assume at the loop base is that P lb s 0. This condition results 
naturally in the chromosphere, irrespective of what is assumed at the base 
as long as the base is plaoed sufficiently deep. If the base is placed 
within several scale heights of the top of the chromosphere, then it is the 
only boundary condition that is appropriate; the use of rigid-wall 
conditions, or any others, is likely to lead to incorreot conclusicns. If 
the base is placed many scale heights deep, then the particular boundary 
condition used becomes irrelevant. 

III. Analysis Results 
a) Perturbation Equations 

In this paper we calculate growth rates and eigenfunctions only for 
the simplified problem discussed by previous authors; we neglect gravity 
and any spatial variation in the loop area or heating function. We discuss 
the effects of gravity and the area variation in a subsequent paper 
(Antiochos and An 1985). For purposes of comparison, we calculate the 
growth rates and modes for several sets of boundary conditions, although as 

we discussed above, the physically appropriate ones are T,, = 0 and 

lb 


Under the special assumptions above, the perturbation equations (2) - 

(6) can be simplified to a form essentially equivalent to the one we used 
previously (Antiochos 1 979 ) s 

1. (y 4_ tJ) + q( x ) M + Xy = n G(x) , <1«> 

where we have defined dimensionless variables; 


x = s/L 


L is the loop half-length. 


y(x) s (T q /T 0 (0)) 7/2 , 


5/14 , 
_ d 


V - 


— V 
V ds 1 


and n = v~ s constant. 

0 

The eigenvalue X is defined as: 


X “ “ (L/L c r v T , 


where x is given by the coronal cooling time scale 


T = f ' >° 6 P 0 ^ (T Q ( 0))" 7/2 


The functions Q(x) and G(x) are given by: 

-*(££) +a ( T o!t) 


G(x) =-y 5/U [ 3/5 X + + a^ Q 9? +aP 0 gp n 


(15) 

(16) 

(17) 

(18) 

(19) 


( 20 ) 


(21a) 


(21b) 


where a is a constant of the dimensions: 
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The boundary conditions oust now be specified. Assuming the "natural" 
ones, ^ = T-y } s 0 leads to: 

n = 0 and = 0 (23) 

at each end of the loop. Note that we have used the fact that the 
equilibrium heat flux vanishes at the base in order to derive 
Equation (23). If one assumes instead the rigid-wall conditions 

V it = T 1k = 0, this leads to the constraint (Antiochos 1979): 

— 5 / 1 A 

p y dx = 0 (24a) 

and the conditions: 

P b + y b 5/14 n = o . (24b) 

With this formulation it is difficult to compare the effects of the two 
sets of boundary conditions since one set leads to homogeneous conditions 
on p, (23), whereas the other does not. In order to facilitate the 
comparison, let us for the moment reformulate the problem in terms of a new 
variable defined as G = p + y 5 ^ 14 n, so that the boundary conditions (24b) 
will become homogeneous. Hence, the problem reduces to solving an equation 
for £ and n of the form (14) (only the form of G(x) changes) and with 
boundary conditions: 

for free surface, £, =0, and r| 

b 

for rigid wall, £, = 0, and n 

b 

The relation between the two sets of conditions is now transparent. We 
note that if, as expected, £ decreases sufficiently rapidly in the 
chromosphere, then (24) implies that in the limit of large chromospheric 
depths, i. e. , L c /L -> 0, 

n -y (L c /L) <s y" 5/l4 > 




(26) 


where the brackets indicate the average ooronal value. Hence, for large 
depths, n decays as L c /L. This is the- same behavior that we found before, 
except that now since there is no gravity, the only scale in the problem is 
the coronal loop length, so that this length rather than the gravitational 
scale height determines the decay scale for the effects of the wall. 


b) Aaifflcfcflfclfi £orm 


In order to check whether 5 will, in fact, decrease rapidly in the 
chromosphere, let us examine the asymptotic forms for ? in this region as 
predicted by equation (111). We first consider the free-surface case, 
Sinoe T q is approximately constant in the chromosphere, Equation (1H) 
reduces to: 

2 2 

>b 72 + 4 o tr) h - (r) T 5 f ■ 0 • < 2 ?> 

dx \ 0/^ ' c ' 

where the subscript B f n is used to refer to the free-surface variables. 
Assuming that the radiative loss function has a simple power-law dependence 
on T, and that the power is positive so that the chromosphere is 
radiatively stable, then we expect that: 



(28) 


where x, is the radiative cooling time at the base. Note that our 


definition for X has a sign change from that used by Field (1965) and that 
the negative sign is required for stability. Assuming that the 
free-surface growth rate, , is of order the coronal cooling rate 1/x, 
Equations (27) and (28) yield the following asymptotic form for c, 

-s/£ 


S ~ e 
a ~ l 


.7M 


- 1/2 


,V°) 


where 


(1 + x/x b ) 


(29) 





The longest length scale for the deoay is obtained In the limit t/t^ -> 0. 
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This length, “ I^T^/Tq^)) * * s tyPi^Hy small compared to physical 

size scales either in the corona or chromosphere. We show below that the 
form predicted above is obtained from exact solution of (14), 

We now oonsider the rigid-wall case. With ? as the dependent 
variable, Equation (1*4) becomes: 


fe {y E ? r) + Q<x> ? r + \ ? r " "r <\ f(x) + S<*» ' 


(30) 


where the subscript n r w is used to indicate rigid-wall quantities. The 
functions f and g are given by: 


f(x) 


, l . 5/14 


(31) 


and 


, . 5/14 / 4 dfy ; 

g(x) =-y I t — k + 
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oP, 


zx 


0 3P, 


(32) 


dx 0^ 

Following the arguments above, in the base region, Equation (30) reduces 

to: „ . . 2 


2 / v c I \ ? 

C( (^) v v { k ! 


T V C 
r r 


T + T/T 


) ' 


(33) 


where, again, due to the form of the radiative losses, we assume that: 



(34) 


The general solution to (33) consists of a linear combination of the 
homogeneous and particular solutions. However, the homogeneous solution is 
equivalent to the free-surface solution, and we found above that this 
decays exponentially. Hence, the solution at large depths is dominated by 


the particular solution: 


Xf, as expected, v 
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2 , 5/M (VVJ«) 

5 n r y b (v T + 1 ) 

t b ' 


0(1/r), then fort/r, -> 0 
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2 „ 5/14 

£r ~ y 5 y b 


(35) 


(36) 


We will see below that (35) also agrees with exact solution of (111). 

IC1 
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These results verify our claim that n decays as L /L. We have that 

L c /L r 1 


W 


_ -5/14 . . 

5 r y dx + 


L c /L 


-5/14 , 
y dx 


(37) 


Substituting i the asymptotic form (35), or (36), into the second integral 
yields (26) to within a factor of order unity. 

The effect on the growth rate of imposing a rigid-wall constraint can 

also be estimated. Using the equations for 5. and S , and their boundary 

f r 

conditions, the following expression is obtained: 



If L^/L « 1, then to lowest order in L c /L, Equation (38) can be 



0 


Each of the integrals in (39) is of order L c /L; hence, 6v is also 
order. 


(38) 


written 


• (39) 

of this 


Depending on the sign of <5v, the rigid-wall can have a stabilizing (6v 
negative) or destabilizing (<5v positive) influence. Since the lowest 


order mode, has no zero crossings and its normalization is irrelevant, the 
sign of 6v is determined primarily by the sign of the terms in the brackets 
in the rightmost integral; however, this factor does not have a unique 
sign. The first term is always positive and for models of interest has 
magnitude ” 1/t; consequently, it enhances instability. The sum of the 
second and third term is proportional to (-e/P^ ) for the case where the 

heating e is constant per particle, and (-2£/Pq) for e constant per unit 

volume. Hence, the contribution from the last two berms is negative, and 
3ince in the static model the energy input and the radiative losses are 

roughly equal in the corona (e.g. , Vesecky, Antioohos, and Underwood 1979), 
g/Pq is also of order 1/t. It is, therefore, not clear which of the two 
contributions will dominate; hcwever, it turns out that the negative 
contribution is usually larger than the positive one so that the rigid wall 
has a net stabilizing effect. 

This is to be expected physically. A rigid wall acts to inhibit flews 

and, consequently, any density disturbance. For optically thin radiation 

at coronal temperatures, both the temperature dependence and the density 

dependence of the losses is such as to promote instability (Field 1965). 

By suppressing the density effect, a rigid wall tends to weaken the 

instability. It is interesting to note that this need not always be the 

case. Sometimes a rigid wall can actually enhance Instability. In 

particular, if is very large and dominates the other terms in (39), we 

would expect 6v to oe positive. This situation does occur for models in 

which the base is placed high up in the transition region and the 

temperature boundary condition is assumed to be dT /ds = 0. For such 

lb 


models the growth time l/v^ is very short, of order the oooling time in the 
transition region rather than that in the corona (Antiochos 1979). We find 
that for these oases the rigid-wall growth rate v r is slightly larger than 
Vf, However, for models in which some cool stable material occurs at the 
base, the rigid-wall boundary conditions always result in smaller growth 
rates, 

iv. numerical RaaaLts 

The mathematical problem that must be solved is now clearly defined. 
For the free-surface conditions, (14) reduces to a standard Sturm-Liouvllle 
problem, as in Antiochos (1979), which we solve in exactly the same manner 
as before. 

For the rigid-wall conditions, the inhomogeneous problem must be 

solved subject to the constraint (24a). Note that this is not a standard 

Sturm- Liouville problem; hence, none of the well-known theorems (e.g. , 
Morse and Feshbach 1953) on the behavior of the elgenfunotions and 
eigenvalues need apply. There are two parameters, X and n , that must be 

determined; therefore, we do not use a "shooting" technique as in the 
homogeneous case to obtain solutions. First, we note that since n will, in 
general, be non-zero for any finite-length loop, we can eliminate it by 
defining a new dependent variable (P/n). Now (14) and (24a) define a 
nonlinear boundary- value problem for (p/n) and the unknown parameter X. 
This is solved by standard techniques: Newton-Raphson iteration and 

finite-difference solution of the linearized equations. Convergence was 
rapid, < 20 iterations for all the cases we investigated. 

One numerical complication is the extreme variation in temperature 
scale height between the corona, the transition region, and the 




chromosphere, This makes it impossible to use a uniform grid for finite 
differencing the independent variable x. However, we cannot use the 
equilibrium temperature as the independent variable since the 
temperature gradient vanishes at the loop apex and bane. Instead, we 
define a now independent variable r: 



and use a uniform grid in r. It is evident from (30) that the effect of 
transforming to r is to increase the number of grid points at transition 
region temperatures (T Q (0) < T Q < T b ) by a factor of order (T fl (0)/T Q ) 2 . 
We found that this was sufficient to yield an accurate solution, 

a) £flnn jq£ Heating and Cooling 

There are two functions in the problem that have yet to be specified: 
the radiative loss coefficient, A ( T) t and the energy input rate e(n,T). 
For the former we use a smooth analytic approximation to the curve derived 
by Raymond, Cox and Smith (1976). Figure 2 shows this ourve along with our 
approximation to it. The stability properties of the models are not 
sensitive to the exaot form of A(T) as long as this form is such that it 
implies radiative instability above " 10 and stability below. 

Since the mechanism for coronal heating is not known, the form of the 
energy input e is essentially arbitrary. In our previous work we found 
that the stability properties are insensitive to the exact form for e, at 
least, for the case where this form is a power-law dependence on n and T. 
Hence, we simply assume that in the corona and transition region the energy 
input either per unit volume or per particle is constant. At lower 
temperatures a different form for the energy input must be used in order to 
have a significant mass of material at "chromospheric” temperature, 


T < 5 x 10^ K. In order for the model to have a large chromosphere, the 
temperature gradients and, hence, the heat flux must be small at these 
temperatures. The radiation losses in the chromosphere cannot be balanced 
by conduction from above (e.g. , Athay 1981); thus for a static model these 
losses must be balanced by the energy input (e.g., Craig, Robb and Rollo 
1982). 

Therefore, we use the following form for the energy inputs 

e = Qn y tanh(x (T)) + n 2 A(1 - tanh(x(T))) , (41 ) 

where 

X (T) = C (T/T b - 1 + <5) m (42) 

and Q,y, C , m and S are constants. Equations (41) and (42) imply that for 

6« is 

1 Qn Y , for T > T fo 

(43) 

n 2 A, for T = T fa , 

v 

i. e. , the energy input rate tends to Qn for temperatures larger than the 
base temperatures and tends to the radiative loss rate at the base 
temperature. The values 0 and 1 for Y correspond to constant coronal 
heating per unit volume and per unit mass, respectively. The depth of 
material at the base temperature is determined by the value of <5 ; it 
becomes arbitrarily large for arbitrarily small 6. 

The form that we assume for the energy input, (41) and (42) may appear 
highly unphysical because it undegoes an abrupt although continuous change 
near T = T b » However, such a sharp structure in the energy input rate is 
implied by the observations. Empirical models of the upper chromosphere, 


ouch as the VAL models (Vernazza, Avrett and Loosser 1981), all Indicate 

a 

the presence of a flat temperature n plateau n at T “ 2 x 10 K with a very 
abrupt rise at higher temperatures. Within the context of the static 
model, and neglecting possible kinetic effects (Shoub 1983), the only way 
to obtain the observed very sharp structure in the temperature profile is 
to have a correspondingly sharp structure in the energy input rate. This 
is one reason why the static models are generally assumed to have a base 
above 2 x 1 0 /f K . The particular form, (41 ) and (42) is chosen simply for 
convenience. It permits us to investigate the effects of adding a varying 
amount of cool, radlatively stable material to the loop base. 

b) firmth .Batea 

(1) Free Surface 

We have calculated the eigenvalues 1 for a wide range of statlo 
models. In particular, we have determined the dependence of the stability 
on the following properties of the model: the coronal temperature, the 

base temperature, the depth of the chromosphere, the form for the heating 
and, of course, the boundary conditions. The main result of this analysis 
is that for the boundary conditions T ib s °* ^ = Of all the models were 
unstable with a growth time approximately equal to the coronal cooling 
time. The growth rate was found to be completely insensitive to the loop 
temperature either in the corona or base, to the amount of chromospheric 
material or to the form of the chromospheric heating. Hence, we conclude 
that the instability is a true physical one and will occur in solar loops 
that satisfy our basic assumptions, i.e. , spatial variations in the loop 
area and heating, and the effects of gravity can be neglected. 

Figure 3 shows the growth rate vt of the lowest free-surface mode 


obtained by solving (14) for a series of statlo models with different base 

* 

properties. These statlo models all have very similar coronal properties. 

The apex temperature T Q (0) was chosen in each case to be 10^ K and the 

-1 2 

coronal heating given by y s 1 and Q = 10 ergs/sec/particle, equation 
(41). Qiven the coronal temperature and heating rate* then as is well 
known (e.g. , Veseoky, Antiochos, and Underwood 1979) * the coronal density 
and coronal loop length are constrained by the scaling laws. They turn out 

q -3 8 

to be Oq(0) s 3.2 x ICr cm and L = 5.2 x 10 cm. Conditions at the 
base, however, are very insensitive to the coronal conditions, so that we 
can vary the base temperature and the amount of base material with almost 
no change in coronal parameters. 

In our model the base temperature is given by T b , and the depth of the 
chromosphere is determined primarily by the constant 6, equation (43). 
Hence, there are actually two distinct sets of models in Figure 3. In one 
set all the models have a value of 6 =. 1 , which implies that there is an 
insignificant amount of material at the base. These models correspond to 
the usual static-loop models, as in Vesecky, Antloohos, and Underwood 
(1979). For such models we investigated the effect of the base temperature 

5 4 

on stability by varying T b from 2x10 K to 10 K. In Figure 3 only the 

4 

cases down to T b s 3 x 10 K are shown, but no change was observed in using 

4 

T b down to 10 . As is clear from the figure, the value of the base 
temperature has essentially no effeot on loop stability. The growth rate 
vx is of order unity irrespective of whether T b is chosen to occur at 
chromospheric or transition-region temperatures. 

The second set of models in Figure 3 are those with a fixed base 

4 

temperature T b = 3 x 10 K, but wi'ih varying amounts of base material at 

—6 

this temperature. The value of 6 was varied from 6 = 1 to 6 =5x10 . 

The other parameters in the energy input rate (equation (42)) were fixed, 


■n 




C b to and m «s 3 , This value for m la convenient numerically since for 
m b 3 the depth of the base region turns out to be proportional to 
1/ 6; hence, the range of chromospheric depths covered by the range UBed for 

^ r 

6 la correspondingly large, (L - L 0 )/L c ranges from " 10 to " 10. Again, 
we note that there is almost no variation in v Irrespective of the depth 
of the chromosphere. 

We have investigated, the effect of the form of the coronal heating on 
stability by performing a series of calculations as in Figure 3 for the 
oase of uniform heating per unit volume, y a 0. The results are 
essentially identical to those shown in Figure 3 . This agrees with our 
previous calculations (Antiochos 1979) in which we investigated a wide 
range of values for Y and found no significant effect on the growth rate. 

We have also looked for any possible effect due to the form of the 
chromospheric heating by varying the exponent m defined in (43). This 
parameter determines the degree of stability of the chromosphere, i.e. , the 
magnitude of the quantity (T o aX/0T Q ) (of,, Field 1965). The chromospheric 
magnitude of this quantity is small in our models compared to its coronal 
value except for the case m s 1 . Note that for m s 1 the chromospheric 

depth L - L c varies only as ln(1/5) so that even with the smallest value of 

“15 

6 usable numerically, "10 , the chromospheric dijpth is not large, 

(L - L c )/L c “ .05. The results for the cases with m s 1 are identical to 
those shown in Figure 3 for the case m = 3 . Again, there is no significant 
dependence of the growth rate v on either the base temperature or the 
chromospheric depth. We conclude that the instability does not depend on 
the particular form, (42), that we assumed for the chromospheric heating. 
This result is to be expected. For temperatures < 5 x 10^ K, the 

quantity — — is strictly negative in all our models; hence, this region 

u ° i o 

oannot contribute to any instability. As we will see below, the unstable 


modes have a vanishingly small amplitude In this region compared to coronal 
amplitudes. 

(11) Rigid Wall 

We have solved Equation (14) for the rigid-wall conditions (24) and 
for the exact wide range of equilibrium models used with the free-surface 
conditions. Some of the results are shown in Figure 3 . The growth rate of 
the lowest rigid-wall mode is plotted for comparison with the corresponding 
free-surface mode. The rigid-wall rate is seen to approach the 

free-surface one as the chromospheric depth becomes large, I^/L << 1 ; the 
ratio of the growth rates is near unity, v r /v,> b ,85 for the deepest model 
in the figure, L c /L = .07. These numerical results confirm the analytic 
arguments above that for large depths the rigid-wall condition is 
equivalent to the free-surface one. 

We find that the rigid-wall growth rates are also insensitive to the 

coronal temperature, the base temperature or the form of the heating. 

However, Figure 3 shows that they do depend on the depth of the base 

region. For models with L c > L/2 , the growth rate is negative, indicating 

stability. Note that this value for !. / L actually implies a large 

chromospheric region. Sinoe the base density is approximately two orders 

of magnitude greater than the coronal density, the bulk of the loop plasma 

is in the base for L c /L = .5. Note also that the damping rate is almost 

•<•2 

constant for models with 6 2 10 . This result emphasizes the unphysical 

nature of the rigid-wall conditions and the danger in using them. If one 
assumes the rigid-wall condition and then considers models with increasing 
chromospheric depths, it first appears that statio loops are stable and 
that the stability is insensitive^ to base conditions. One would naturally 
conclude that the models are physically stable. But as the base depth is 


increased, the models suddenly become unstable. Somehow the addition of 
too much stable material at the base destabilizes the coronal This is 
clearly an unphysical result and is due to using improper assumptions at 
the onset. 

The same type of behavior is found when boundary conditions other than 
T lb = 0 and n = 0 are used. For example, Craig and McClymont (1981) find 
that the conditions dT^/ds s 0 and P s 0 lead to abrupt changes in the 
growth rate at very small chromospheric depths. We have also examined 
these conditions and a variety of others defined by the vanishing of a 

linear combination of T and F or P„, and V . Depending on the 

lb lb lb lb 

particular equilibrium model, the boundary conditions assumed can result in 
growth rates either larger or smaller than the free surface one; however, 
in all cases the rates tended to the free- surface result for large 
ohromospheric depths. In all but the free-surfaee case the growth rates 
exhibited an abrupt change at some range of chromospheric depth; hence, the 

conditions T s n =0 are the only ones that are generally valid. 

lb lb 


c) Normal Modes 


Along with the eigenvalues 

A we have 

also 

calculated 

the 

eigenfunctions p for 

all the 

equilibrium models 

considered. 

TWo 

representative oases are 

shown in Figures 4a and 4b. 

In 

Figure 4 a we 

show 

both the free-surface 

Tj/Tq and 

the rigid-wall 

mode, T^Tq, for 

the 


equilibrium model of Figure 3 with the deepest chromosphere (L /L = .07 

c 

corresponding to 6 = 5 x 10" 6 ). In Figure 4 b the results for the model 
with 6 s 1 0 ” 3 , (L c /L s .97), are shown. 

We first discuss the results for the deepest model, Fig. 4a. In order 
to resolve their structure in the transition region, the modes are plotted 
against a nonlinear distance scale. From the loop apex, s = 0, to a point 






two coronal depths down, s = 2L C , the equilibrium temperature T Q is used as 
the independent variable. From this depth down to the loop base, s = L, 
the distance itself is used. Of course, all the significant structure in 
the eigenmodes occurs well above the point s = ZL^. The region below this 
depth is included in the figure only to point out that the numerical 
calculations verify the analytic asymptotic forms obtained above. Using 
(29) and taking t/x b s 0, we find that the decay scale A for the 

_3 

free-surface mode is equal to 2.16 x 10 L c> Hence, the free-surfaoe mode 
should have essentially zero amplitude in the deep region. The numerical 
results agree with this. We find that the amplitude of T^ /T Q drops below 
computer accuracy well before s = 2L . On the other hand, the amplitude of 
the rigid-wall mode, T^ /tiT q , does not vanish. Instead, it has a constant 
value of 0.4 over most of the chromosphere and vanishes only near the base 
on the size scale JL. This behavior is also in agreement with the analytic 
results. In the chromosphere, T x /n T Q s 5 /n y^ 1 *; therefore, from ( 36 ) 

1 

we expect the mode to approach the value 2/5, except very near the base, 
where it must vanish due to the boundary conditions. 

It can be seen from Fig. 4a that the modes have significant amplitude 
only in the corona and transition region. The free-surface mode has been 
normalized to equal the rigid-wall value at the apex. Note that the 
rigid-wall quantity T^ /r\ T Q does not have an aribtrary normalization 
factor since this quantity is the ratio of the temperature- to- pressure 
perturbations. Its magnitude has physical significance. Our results 
indicate that the pressure perturbation is in the same sense but smaller 
than the temperature perturbation. As the depth of the chromospheric 
region increases, L c /L -> 0, we find that the coronal magnitude of T^/hTq 
increases as (L c /L) . This confirms our previous result that n should 
decrease as L c /L for a fixed temperature perturbation. 


On the scale of Pig, 4a the two modes are almost identical In the 
corona and transition region. Only near the base is there a significant 
deviation; and as the depth of the base region is increased, the two modes 
become indistinguishable. The modes exhibit a stong peak near 5 x 10^ K. 
At approximately this temperature the radiative losses change from being a 
destabilizing to a stabilizing effect, so that below this temperature the 
modes decay very rapidly on the scale of It is important to emphasize 
that although the modes peak strongly in the transition region, the 
instability is not a local one. The modes have finite amplitude throughout 
the corona, and the growth rates correspond to the coronal coaling time 
rather than the much faster cooling ( in the transition region. Hence, 
the instability is a global one involving the complete loop except for the 
chromosphere. 

The situation is somewhat different for the model with a shallower 
base, Fig. 4b. The free-surface mode is very similar to that of the model, 
with the deeper base. Indeed, this mode is insensitive to just about 
everything; hence, the curve shown in Fig. 4a (or 4b) can be considered the 
universal form. 

However, the rigid-wall mode is obviously very different from that of 
Fig. 4a. The most striking difference is that the mode has a zero crossing 
at the base of the transition region and is negative in the chromosphere. 
This is in agreement with the analytic results. We note from (35) that as 
v decreases from a large positive quantity, the asymptotic amplitude 
becomes negative. As v decreases further to a large negative quantity, 
the amplitude becomes positive again. The numerical results agree with the 
behavior. For models with large depths, 6 < 10"^ , the amplitude in the 
chromosphere is positive. For models with intermediate depths, 6 “ TO -3 , 
it is negative. It becomes positive again for small depths, 6 > 10“^. For 


models without a very deep base, the chromospheric form of the rigid-wall 
mode is clearly quite sensitive to the exact position of the wall. The 

free-surface mode, on the other hand, shews no ohange for all values of 6 

that we Investigated. 

In the corona the rigid-wall mode has the same general shape as that 
in Pig. 4a, but the amplitude is an order of magnitude smaller. This means 
that for a given amplitude temperature perturbation, the effect of the 
decrease in the base depth is an increase in the coronal pressure 

perturbation by an order of magnitude. In fact, the pressure perturbation 

actually has a larger coronal amplitude than the temperature perturbation. 
We note from Fig. 4b that the ratio of the temperature to pressure 
amplitude is less than unity for temperatures down to ' 8 * 10 K, This 
region encompasses about 70 percent of the coronal portion of the loop. 
Since from (5) 


T 1 n i 

v + v >‘ 1 


(W 


our results imply that over most of the oroona the temperature and density 
perturbations must be in the same sense. A decrease (an increase) in the 
temperature is accompanied by a decrease (an increase) in the density as 
well. Since a decrease in the temperature implies an increase in the 
radiative loss rate (for T > 10 5 K), whereas a decrease in the density 
implies a decrease in the losses, the perturbations oppose each other as 
far as instability is concerned. The net result is stability, v <0 for 
the model in Fig. 4b. For models with a deeper base region, as in Fig. 4a, 
the temperature perturbation is much larger than the pressure perturbation, 
and, consequently, the density perturbation takes the opposite sense. Both 
the temperature and density disturbances act to promote instability, and 
the model of Fig. 4a is, indeed, unstable, v > 0. We conclude that the 


stability of the rigid-wall models such as Fig. Mb is a boundary-condition 
effect. It is due solely to the assumption of a rigid-wall at a particular 
depth. 

V. Elfiflttafllgn 


There are 

two key 

results of this paper. 

The 

first is 

that 

the 

appropriate 

boundary 

conditions to assume 

at 

the loop 

base 

are 


T ib = n lb = 0. For models with a sufficiently deep base 

(L, » min (H.,L )), all boundary conditions can be usedj however, the set 

b be 

above is the only one that is appropriate for models that do not satisfy 
this condition. Indeed, we find it surprising Just how universally valid 
these boundary conditions appear to be. Even for models with little or H2 
chromosphere, they yield the correct growth rate and form of the eigenmode 
for the instability. Contrary to the conclusions of several authors 
(Habbal and Roaner 1979; Chiuderi, Einaudi, and Torricelll-Ciamponi 
1981; Craig and McClymont 1981; and McClymont and Craig 1981a, b, c), the 
mere presence of radiatively stable material at the loop base does not 
stabilize the static models. If incorrect boundary conditions are used, 
then the base can act either to stabilize or, as we showed above, to 
destabilize the static models. However, if the proper conditions are used, 
then the amount of base material has no effect on the stability. 

This is a convenient result from a numerical point of view. It is 
often quite difficult to handle numerical models with extremely deep base 
regions (e.g. , Craig, Robb, and Rollo 1982). Our results indicate that 
models with very shallow bases can yield correct results if the correct 
boundary conditions are used. 

The second main result of this paper is that under our assumptions, 



static models of coronal loops are thermally unstable. The growth rate of 
the instabllty is of the order the coronal cooling time (20), The basic 
assumptions that we made for the equllbrium models were that gravity and 
all spatial variations in the loop area or heating rate are negligible. 
These are the most commonly used assumptions in calculating static models. 
An important question is: how sensitive are our results to these 
assumptions? MoClymont and Craig (198la,b, o) have considered the question 
of spatial variations in the heating rate. They find that these can act 
either to enhance or to damp instability. Lacking a model for the coronal 
heating process, nothing definitive can be stated on this issue; hence, we 
will not consider it further. We will discuss the effects of gravity and 
the area variation in a forthcoming appear. In general, the area variation 
acts to enhance the Instability; however, we do not expect it to enhance 
the growth rates to values much beyond the ones here, since these are 
already of the order of the coronal cooling rate. Gravity, on the other 
hand, acts to damp the instability (e.g. , Wragg and Priest 1982); hence, 
our results are certain to apply only for loops in which gravity can be 
neglected. 

We can Identify at least two possible physical situations in the solar 
corona where the effects of gravity should be small so that our results 
clearly apply. One is the case of loops sufficiently low lying, heights 
K 1000 km, so that the loop height is smaller than any gravitational scale 
height in the loop. Another is the case of magnetic field lines that are 
concave in the corona, so that near the loop center the force of gravity is 
directed toward the center rather than away from it, as in the usual case. 
This is the type of magnetic geometry believed to be responsible for 
quiescent prominences. Our results imply that these two types of 
structures are naturally unstable without any special requirements on the 


coronal heating process. We believe that the thermal instability discussed 
in this paper is responsible for the cool material that is often observed 
to occur in the corona. 

Of course, in order to prove this statement, the nonlinear development 
of instability must be calculated. Clearly, the instability must not 
saturate at a low amplitude if it is to produce observational effects in 
the corona. Only the linear growth rates were considered in this 
paper; the nonlinear evaluation must be determined by numerical simulation. 
Numerical simulations of coronal loops have been performed by several 
authors (e.g. , Craig, Robb, and Rollo 1982; Peres .efc al. 1982; Oran, 
Mariska, and Boris 1982). These simulations have not exhibited any 
significant evidence for thermal instability; however, to our knowledge no 
one has considered models in which the effects of gravity are negligible 
and the correct boundary conditions are employed at the loop base. Our 
results indicate that the instability should be present in such models. We 
intend to investigate the nonlinear stability of coronal loops in a future 
work. 
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Flftur.fl Captions 

Figure 1 Plots of free-surface and rigid-wall pressure perturbation 
(Figure la) and mass flux perturbation (Figure 1b) for the 
highly simplified model described in the text. The 
free-surface results are Indicated by solid ourves and the 
rigid-wall one by broken curves; the two types of modes are 
distinguishable only near the base, s = L. 

Figure 2 The radiative loss coefficient of Raymond, Cox and Smith (1976) 

(broken line) and our analytic fit to it (solid line). The units 
in the figure are arbitrary. 

Figure 3 The free-surface (solid curve) and rigid-wall (broken curve) 
growth rates for a series of statio models with varying base 
temperatures and/or varying base depths. 

Figure 4 a Plot of lowest free-surface temperature perturbation (^/Tq) 

and rigid-wall perturbation (^Pg/TgP^) for the equilibrium model 
with 6 s 5 x 10 ^ and a length L = 14.2 L^. The free-surface 
mode is indicated by a solid curve, and the rigid-wall one by a 
broken curve; they are indistinguishable in the corona and 
transition region. The scale of the abscissa is log (T q ) for 
0 i s i 2L , and s/L for 2L < s £ L. 

C G 

Figure 4b Same as in Figure 4a but for the equlibrium model with 6 = 10 " 3 , 
and with an abscissa of log (T Q ) for 0 i s/L i .94 and 
log (s/L) for .94 £ s/L < 1. Note that the rigid-wall 
amplitude is negative in the region s/L > .94. 
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